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Abstract 

When the last electron-photon scattering takes place in a magnetized environment, the 
degree of circular polarization of the outgoing radiation depends upon the magnetic field 
strength. After deriving the scattering matrix of the process, the generalized radiative 
transfer equations are deduced in the presence of the relativistic fluctuations of the geometry 
and for all the four brightness perturbations. The new system of equations is solved under the 
assumption that the incident radiation is not polarized. The induced V-mode polarization 
is analyzed both analytically and numerically. The corresponding angular power spectra are 
calculated and compared with the measured (or purported) values of the linear polarizations 
(i.e. E-mode and B-mode) as they arise in the concordance model and in its neighboring 
extensions. Possible connections between the V-mode polarization of the Cosmic Microwave 
background and the topological properties of the magnetic flux lines prior to equality are 
outlined and briefly explored in analogy with the physics of magnetized sun spots. 
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1 Motivations and goals 



The circular polarization of the Cosmic Microwave Background (CMB in what follows) is not 
the direct target of forthcoming experimental searches. It will be argued hereunder that more 
accurate spectropolarimetric measurements of the V-mode polarization can be enlightening 
especially as a diagnostic of the magnetization of the pre-decoupling plasma. The primary 
goal of experimental endeavors in the near future is related, in one way or in the other, to 
the determination of the angular power spectra of the intensity and of the linear polarization 
of the CMB radiation field. Even the B-mode polarization, one of the primary objectives of 
diverse experimental programs, will be unable to shed light on the circular polarizations of 
the CMB. To avoid possible misunderstandings on this point it is desirable to introduce the 
relevant conventions on the Stokes parameters of the radiation field 

/= \E-e 1 \ 2 + \E-e 2 \ 2 , V = 2Im[(£ • e x )*{E ■ e 2 )], (1.1) 
Q=\E-e 1 \ 2 -\E-e 2 \ 2 , U = 2 Re{{E ■ erf (3 ■ e 2 )}, (1.2) 

where e.\ and e 2 are two mutually orthogonal directions both perpendicular to the direction of 
propagation n which is chosen to lie along £3. The temperature autocorrelation (i.e. the TT 
angular power spectrurrjfl) stems directly from the brightness perturbation of the intensity 
of the radiation field: the fluctuations of the space-time curvature act as sources of inho- 
mogeneity of the intensity. The treatment of the polarization is slightly more cumbersome 
and it has to do with the transformation properties of the Stokes parameters. Consider a 
rotation of the unit vectors e\ and e 2 in the plane orthogonal to h and suppose that the 
rotation angle is (p. By rotating the axes in the right-handed sense, i.e. as 

= cos<£>ei + sin<£>e 2 , — — sin<£>ei + cos <pe 2 , (1.3) 

the Stokes parameters of Eq. (11.11) are left invariant (i.e. I' = I and V = V) while Q and 
U (introduced in Eq. (11. 2p ) transform as 

Q' = cos 2cpQ + sin 2(p U, U' = - sin 2<p Q + cos 2(p U. (1.4) 

In different words, while / and V are invariant under a two-dimensional rotation, Q and U do 
transform and do mix under the same rotation. It is worth stressing that he transformation 
parametrized by Eq. (11. 3p is not a global rotation of the coordinate system but it is rather 
a rotation on the tangent plane to the two-sphere at a given point. The transformation 
properties of the Stokes parameters under Eq. (jl.3p allow determine their associated spin 
weight [U El El U\- The brightness perturbations associated with I and V (i.e. Aj and 

2 From the definitions of the Stokes parameters it follows (see, e. g. [U El 13]) that I 2 > Q 2 + U 2 + V 2 , 
where the equality sign arises in the case of the field of a plane wave. 

3 The specific definitions of the relevant angular power spectra are provided hereunder, see, e.g. Eqs. 
(fTTTjl and (Ti~T2|) 
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A v ) have both spin-weight 0. The brightness perturbations of Q and U (i.e. Aq and Ay) 
transform as a function of spin- weight ±2, since, from Eq. (jl.4p . A±(n, r) = Aq (n, t) ± 
«Au(n, r) transform as, A±(n, r) = e T2iI,p A-t (n, r). Consequently, while A±(n, t) have to be 
expanded, on the two-sphere, in terms of spin ±2 spherical harmonics ±2Ye m (n) 

A±(n,r) = 5Z a ±2,£m±25 / £m(n), (1.5) 

Ai(n, r) and Ay(n, r) have to be expanded on the basis of spherical harmonics Y( m (h) as: 
Ai(n, r) = ^ aJJ, ^ m (n), A v (n, r) = ^ 4m ^m(n). (1.6) 



Both spin-0 and spin- ±2 spherical harmonics arise naturally as Wigner matrix elements [5] 
depending in general upon three different eigenvalues: while the third eigenvalue is for 
scalar harmonics, it is =Fs for spin-s weighted harmonics. The complete information on the 
radiation field of the CMB should therefore stem from the analysis of the TT and VV angular 
power spectra 

( TT) _ 1 (n<n*„® \ r (VV) - 1 V (n {V) *n {y) \ M 7) 

1 ~ 2£ + 1 ' a f™,n W ~ 2£ + 1 ' f-mh K 1 -') 

m=—£ m=—£ 

as well as from the angular power spectra of the E-mode and B-mode autocorrelations 0. 
The E-mode and B-mode autocorrelations are defined as 

^(ee) _ 1 /„( E )*„( E )\ r^ BB ) - 1 V" /„( B )*„( B )\ (^ z\ 

~ 2£+l a imh W — 2f _L I \ a lm a lmh \ L -°) 

m=—£ m=—£ 

where and linear combination of the coefficients a±2 £ m already introduced in 

Eq. (USD: 

( E ) _ 1 I i \ (B) « / s /-, n x 

ffm — — -{0-2, £m + «-2,£mJ, <fy m — - (^2, £m ~ Cb-2, im ) • U- y J 

In real space the fluctuations constructed from and have the property of being 
invariant under rotations on a plane orthogonal to n and, up to an ^-dependent prefactor, 
they can be expanded in terms of (ordinary) spherical harmonics: 

A E (n, r) = £ A7 1 a£ Y, m (n), A B (n, r) = £ A7 1 a£> Y tm (jh), (1.10) 



where iV^ = y (£ — 2)\/(£ + 2)!. Before matter radiation equality the radiation field is cus- 
tomarily assumed to be unpolarized. The properties of electron-photon (and ion-photon) 



4 This statement on the autocorrelations of temperature and polarization does not exclude the possibility 
of discussing and measuring the various cross-correlations between temperature and polarization. 
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scattering imply that the radiation can become linearly polarized provided the incident 
brightness perturbations have a non-vanishing quadrupole moment, i.e. A12 7^ 0. 

Circular dichroism as well as linear polarization of the CMB becomes theoretically plau- 
sible in the presence of pre-equality magnetic fields [9j [10] but, so far, there has not been 
any specific calculation of the circular polarization induced by a magnetized plasma prior to 
recombination and in the framework of the concordance model. One of the purposes of the 
present article is to fill such a gap. While the V-mode polarization is suppressed in compar- 
ison with the E-mode polarization, the question is to determine quantitatively the nature of 
the suppression and its typical range in multipole space. The answer to the latter question 
can only be dynamical: it will be interesting, for the present purposes, to understand how, 
when and to what extent a radiation field which is originally unpolarized prior to equality 
will become circularly polarized after photon decoupling. 

Some of the phenomenological aspects of the present considerations can be understood 
in analogy with solar spectropolarimetry. The knowledge of large-scale magnetism is often 
inspired (and partially modeled) by our improved understanding of solar magnetism (see, 
for instance, |TT] for a dedicated review on the subject). Circular polarization is known to 
occur in sunspots [121 DSl EE] an d it is sometimes argued that the analysis of polarization 
might improve the understanding of the local topology of magnetic flux lines in the vicinity 
of sunspots pTj. The amount of circular polarization depends, in the case of the sun, 
both on the gradients of the velocity field and as well as upon the magnetic field topology 
and intensity. A naive approach to the problem would suggest that, absent any velocity 
gradient, V ~ 9 x 10~ 12 -B spo t A / where A is the wavelength in units of nm (lnm = 10 _9 m) 
and -Bgpot is the magnetic field in units of gauss. As noted long ago (see e. g. [121 f!3] ) 
for typical values of the magnetic field (i.e. -B S pot ~ kG) and for optical wavelengths (i.e. 
380 nm < A < 750 nm) we would get V/I ~ 10~ 5 which is one or two orders of magnitude 
smaller than what approximately observed. The latter discrepancy suggested (see e.g. [15] ) 
the important role played by velocity gradients whose contribution should be included for a 
consistent interpretation of the observational data. 

Circular polarization arises naturally also in synchrotron emission [121 dZl [TS] (see also 
[T9"l [2D]). In this second context the amount of circular polarization is generated from the 
linear polarization by the so called Faraday conversion effect. In the case of the CMB the 
relativistic effect associated with the synchrotron emission is difficult to realize. On the 
contrary the pre-decoupling plasma is rather cold and the electrons are non relativistic. 
Furthermore, it is difficult to justify that the initial radiation field should be linearly polar- 
ized already well before equality. The only way Faraday conversion could lead to circular 
polarization is as a secondary effect when the linearly polarized CMB photons impinge on 
the relativistic electrons in the cluster magnetic field |2T] . In this paper, on the contrary, 
the target is to compute the amount of circular polarization induced at early times by a 

5 As usual Aw denotes the £-th multipole of the intensity of the radiation held. With the same notation 
we will refer, when needed, to the multipoles of other observables. 
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pre-decoupling magnetic field and in terms of unpolarized initial conditions of the Einstein- 
Boltzmann hierarchy. To introduce some quantitative (but still general) considerations it is 
practical to define, for immediate convenience, 

^(TT) = 1(1 + 1) C (TT) ^(VV) = £(£ + l) c {W) ^ n , 

2n 2n 

(EE) 1(1 + 1 ) r y(EE) ^(BB) _ £{£+ l) ^(BB) {-\-\o\ 
Lr e — Ly £ , Lx t — U £ , l 1 - 1 ^ 

measuring, for a given observable, the angular power per logarithmic interval of i. With 
the same notation the cross-correlations between different observables can be defined, for 
instance, as G^ E \ G^ VT ' ) and so on and so forth. In terms of Eqs. (11. lip and (11. 12ft the 



typical orders of magnitude angular power spectra can then be be summarized as 

Gf T) ~ 0(5 x 10 3 )(/xK) 2 , Gf E) ~ 0(150) (/^K) 2 , Gf E) ~ C(50)(^K) 2 . (1.13) 

The TT and TE correlations have been accurately assessed by the WMAP collaboration 
[221 I23T |2l] (see also [251 126]). Interesting measurements on the EE correlations have been 
reported, for instance, by the QUAD experiment [291 EQl EH [32] . Other measurements on 
the TT correlation for large multipoles (i.e. I > 1000) have been reported by the ACBAR 
experiment [271 EE] ■ While all the current experimental data are consistent with the standard 
ACDM paradigm they can also be used to estimate the parameters of a putative magnetized 
background. In this respect, the result is that large-scale magnetic fields of nG strength 
and slightly blue spectral indices are allowed by current CMB data [331 El]- The estimation 
of the parameters of a magnetized background led to the first estimate of the likelihood 
contours in the two-dimensional plane characterized by the magnetic spectral index and by 
the magnetic field intensity. In a frequentist perspective the results of [331 El] exclude, to 
95% confidence level, a sizable portion of the parameter space of magnetized models centered 
around comoving field intensities of 3.2 nG and magnetic spectral indices^ n-Q = 1.6. The 
latter results hold when the underlying model is just the ACDM (concordance) paradigrrjf] 
but the addition of, for instance, dark energy fluctuations does not change quantitatively 
the exclusion plots [IT] . 

The B-mode polarization has not been measured yet and upper limits exist by various 
experiments. In the standard ACDM paradigm with no tensors the BB power spectrum 
vanishes. For instance magnetic fields of (9(5 nG) and blue spectral index lead, at inter- 
mediate scales, to an angular power spectrum which can be as large as 10~ 3 (/iK) 2 |42j. 
The latter estimate can be compared, for instance, with the BB angular power spectrum 



6 For sake of simplicity the numerical values quoted here refer to the maximum of each power spectrum 
7 The conventions on the magnetic spectral indices n-Q are exactly the same as the ones employed for the 

spectral index of curvature perturbations (denoted as n s ). The scale-invariant value corresponds, in both 

cases, to 1. 

8 In the acronym A stands for the (non-fluctuating) dark energy component while CDM stands for cold 
dark matter. 
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expected from the tensor modes of the geometry and from the gravitational lensing of the 
CMB anisotropies. The tensor modes are the conventional (potential) source of B-mode 
polarization in the simplest extension of the ACDM paradigm. Defining rx as the tensor to 
scalar ratio, a typical value rx ~ 0.3 would imply at intermediate multipoles (i.e. £ < 100) 
a BB angular power spectrum of the order of 10~ 2 (/uK) 2 . 

The degree of circular polarization computed here depends, in the simplest case, upon 
the amount of curvature perturbations, upon the magnetic field parameters and upon the 
typical frequency of the experiment. For nG magnetic fields and for a reference frequency 
O(10) GHz the VV angular power spectrum is G^ V ^ ~ 10~ 15 (/iK) 2 . For the same range 
of parameters the cross-correlation G^ VT ' ) ~ 10~ 6 ' 5 (/iK) 2 . If the observational frequency 
decreases the signals can be larger. While this comparison will be more carefully performed 
in section 4 it is important to appreciate that, in some sense, the absolute magnitude of 
the different correlation functions represents just a necessary but insufficient guide for the 
observer since the systematics associated with the circular polarization are different from 
the ones arising in the case of a linearly polarized signal [35J. The second point we wish to 
stress is that, as it will be apparent from the subsequent analysis, low frequency instruments 
seem to be experimentally preferable [361 EZ] an d, in this respect, it is tempting to speculate 
that very low frequency (radio) techniques could be appropriately adapted [381 E3 HQ] ■ The 
third point related to the potential observations of the effects is that spectropolarimetric 
techniques should be probably employed given the necessity of a simultaneous determination 
of the brightness perturbation of the intensity and of the V-mode polarization (see, in this 
respect, section H]). 

The analogy of the present problem with the physics of the sunspots suggests the possi- 
bility of connecting the amount of (measured) circular polarization with the topology of the 
magnetic flux lines. This topic is, in principle rather rich and, in this paper, we will merely 
scratch the surface by presenting some particular examples which are only semi-realistic and 
which are borrowed from known examples in plasma physics. Indeed, the study of the topol- 
ogy of magnetic flux lines has a long history going back to the pioneering work of Fermi 
and Chandrasekhar on the gravitational stability of the galactic arm [4"3"| HI] . It would be 
interesting, in perspective, to connect the possible occurrence (or absence) of the circular 
polarization with magnetized plasmas which minimize the magnetic energy while the helicity 
is conserved (as it should) at high conductivity (see, along this line, the seminal papers of 
Chandrasekhar, Kendall and Woltjer |4"5 ] 136 1 147] ) . 

The layout of the paper is the following. In section [2] the photon electron and photon 
ion scattering will be discussed in the presence of a magnetic field in the guiding centre 
approximation. Details are also reported in appendix |X] In section [3] the same problem will 
be addressed in the case of a magnetic knot, i.e. a simple example of static configuration 
minimizing the energy at a fixed value of the magnetic helicity. It will be speculated that the 
degree of circular polarization could be eventually connected with the topological properties 
of the magnetic flux lines in the plasma. In section H]the evolution equations of the brightness 
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perturbations will be deduced and solved both analytically and numerically. Relevant details 
on this topic are given in appendix [B] The tight-coupling approximation will be applied to 
the new framework and analytical results for the V-mode power spectra for large angular 
scales will be derived. For smaller angular scales numerical results will also be presented and 
compared with the temperature and (linear) polarization anisotropies. Section contains 
the concluding considerations. 



2 Magnetized electron-photon scattering 

The electron-photon scattering is customarily computed without taking into account the 
contribution of the magnetic field itself to the scattering matrix. This happens not only in 
the case when magnetic fields are assumed to be absent but also in the presence of large-scale 
magnetic fields (see, for instance^, [63l|6lJ[65]). The purpose of this section is to drop such an 
assumption and to derive the appropriate scattering matrix for electron-photon scattering in 
a weakly magnetized plasma. It is practical to define, for the present purposes, the outgoing 
and ingoing Stokes vectors whose components are the Stokes parameters, i.e. 



Z out ( W , n, //, <p, tp') = (h, h, U, V), (2.1) 

r>, //, <p') = 4 u', V), (2.2) 

jj, = cos "& and y! = cos$'. The intensity / and one of the components of the linear po- 
larization (i.e. Q) have been replaced, as usual, by I\ = (I + Q)/2 and I 2 = (I — Q)/2. 



The components of the ingoing Stokes parameters have been distinguished by a prime and 
they depend upon // = cos^' and tp'. The Stokes parameters depend upon the (angular) 
frequency u = 2nv. By definition, the scattering matrix connects the outgoing to the ingoing 
Stokes parameters as: 

ZT*(w, M> /A <P, V 5 ') = Sij(u, fi, fi', <p, p')T^(uj } fi', (f'). (2.3) 



The coordinate system has been fixed a; 




f = (cos (p sin i?, sin (p sin i?, cost?), (2.4) 
9 = (cos p cos sin^cos'i?, — sin?9), (2.5) 
p = (— sin <p>, cos ip, 0). (2.6) 

9 It is appropriate to stress that it would be rather pretentious to give complete and thorough list of 
references in connection with primordial magnetism. The easiest solution is to refer the interested reader to 
the dedicated review articles of [9l [10] where a more complete bibliography can be found. 

10 The orientation of the coordinate system corresponding to Eqs. (|2.4[) - (|2.6[) implies that i? X ip = f. In 
other classic references such as [2] the orientation is such that d x tp = —f; this different choice entails a 
modification of Eq. (I2.6[) . i.e. tp = (sintp, — costp, 0). The conventions spelled out by Eqs. (|2.4[) - (|2.6|) will 
be followed throughout the paper. 
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The purpose is to obtain the scattering matrix of electron-photon scattering (see, e.g. [2]) 
but in the presence of a magnetic field and in the guiding centre approximation [H] which is, 
in practice, a controlled expansion in gradients of the magnetic field intensity. The derivation 
of the various components of the scattering matrix is reported in appendix |A] In what follows 
only the results will be reported and discussed. Defining as r e = e 2 /m e the classical radius 
of the electron, the various components of the scattering matrix can be written as: 



S 



11 



<Sl2 
<Sl3 

Oi4 
$21 

S22 
<!>23 
S2A 
031 

£33 



e 

2r 2 



2A 3 M(l-^)(l-// 2 ) + C 2 M/iV 



.2 ../2 



4CMA!(u;)A3(u;)/i// yjl - /x 2 ^1 - /i' 2 cos {iff - ip) 



+ C>) 



Al(co) - /»A 2 (u;) 



ji 2 '/j,' 2 cos [2((f' — <f)] 



H 2 C(u)r, 



2r 



Al(u) + f 2 (u)A 2 2 (co) 



r e 2 / e HCHA 2 (a;) 



(2.7) 

c-os[2(/-^,(2.8) 
(2.9) 



A 2 (u,)-/»A 2 (u,i 
C(w)A 1 (a;)A3(a;)/i V /l - /x 2 yfl - /x' 2 sin (if' - ip) 
A 2 H-/ 2 (a;)A 2 M] sin [2(y/ -?)]}, 

H C(w)A 1 (u;)/i / u , - A 3 (u)yjl-fi 2 y/l-fi' 2 cos (^ - (2.10) 
A 2 (u;) + />)A 2 (u,) 

cos [2(^' - v?)] \, 



C»// 2 r 



V J2 
e 



2r 2 

.2/, A r2 



A?( W ) - /»A 2 M 
C 2 Mr 2 



2r 2 
C 2 Mr 2 

2r 2 



A 2 (u;) + /»A 2 (u;) + 



A 2 ( W ) - /»AiM 



A 2 (c) - / 2 (c)A 2 (a;) 



//'sin [2(<p' - if)], 



-f/.(w) C 2 H AiCw) A 2 (a;y, 



^CMAiMAgM^/l-z/V Vl-^sin - 



+ C 2 M//// 



/2 



A 2 (u,) - / 2 (u,)A 2 (u,) 

a 2 m-/ 2 ma 2 m 



sin [2(f'-f)] 
\i sin [2(y/ — if)] 



|-A x (a;)A3(a;)yi -fi 2 y/l -/i' 2 cos (^' - y>) 



+ CH 



a 2 m-/ 2 ha 2 h 



////'cos [2(<p' - If)] \, 



'31 



-^|/ e (a;)C(a;)A 2 (a;)A3(a;)^/l - //^l - //' 2 sin (</?' - ip), 



(2.11) 

cos [2( J, (2.12) 
(2.13) 
(2.14) 

(2.15) 
(2.16) 

(2.17) 
(2.18) 
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041 
^44 



C(w)Ai(w)////' 



-A3(w)^l - /i 2 Vl - /i' 2 cos - \, 
r 2 e f e (oo)C 2 (co)A 2 (uj)f, 



^f e (u)((u)A 2 (u)A 3 (u)Jl - /iVl " M ,2 sin(^ - y>) 



^WC 2 (^> 



A?( W )+/?( W )A*(a;) 



C(a;)Ai(w)A 3 (u;)^/l - /^V 1 " ^ cos " ¥>)}■ 



(2.19) 
(2.20) 
(2.21) 

(2.22) 



In Eqs. f l2.7p -f l2T22|) various (frequency dependent) quantities have been introduced, namely 
fe{oj)-> Ax (a;), A 2 (w) and A3 (a;). Their explicit expressions are: 



Ai(w) 
A 2 (w) 
Aa(w) 

CM 



1 + 



p 1 



P e 
2 



^ 2 ~ g§ c 
W 2 - W| ; 



P 1 



1 + 



P e 
2 

P i 



UJB_ 



- <4 e 



P c 



/.(") 



(2.23) 
(2.24) 
(2.25) 
(2.26) 



where wb e , i an d w p e , i are the Larmor and plasma frequencies for electrons and ions, namely: 



eB 



m P a 



wb i 



eB 



m p a 



CO 



P e 



47r e 2 n 
V m e a(r) : 



p 1 



\ 



47r e 2 n 



(2.27) 
(2.28) 



where B is the modulus of the magnetic field intensity coinciding, in practice, with the lowest 
order result of the guiding centre approximation. As discussed in appendix |Aj Eqs. ( I2.27P 
and ( I2.28P take into account the redshift dependence of the frequency. The magnetic field 
appearing in Eq. ( I2.27P is the comoving magnetic field (see appendix and this means, in 
particular, that the relation to the physical frequencies is given by 



eB 



a(r) = wjff a(r), 



(2.29) 



and similarly for the other quantities of Eqs. (I2.27P and (12.28)) . The logic behind the 
functions defined in Eqs. f!2.23p . (12.241) and (12.25)) is that we want to factorize the electron 
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contribution by keeping track of the ions. According to this strategy, to leading order in the 
ion contributions, Ai(a>), A 2 (w) and A3 (a;) turn out to be: 



A 2 (w) 
A 2 (w) 



\mzy 



l + O 



c 
c 



l + O 



(2.30) 



In the limit / e (w) — > and by correspondingly setting A 1(0;) = A 2 (w) = A 3 (w) = 1, Eqs. 
(I'2.7jl — (I2.22p reproduce exactly the results of Ref. [2] modulo the different orientation of 
the coordinate system (see Eqs. f)2.4l) - fl2.6p and footnote therein) which leads to an overall 
sign difference in the matrix elements containing the sines. In the limit A\(u) = A2(oj) = 
A 3 (a;) —> 1 the present results coincide with Ref. jH] (see also [50] modulo the typo pointed 
out in [1H])- The evolution equations for the brightness perturbations can be written, in 
general terms, as 



+ e'Ai + n%Ai = ip' - n%<p + e'fxv b + dCi(ji, u), 
A Q + e'A Q + n^Aq = e'C Q (/i, w), 
A^ + e'Au + n^Au = e'C v {^ u), 
A' Y + e'A v + rfdiAy = e'CVGu, u). 



(2.31) 
(2.32) 
(2.33) 
(2.34) 



Concerning Eqs. fl2.3ip - fl2.34p few comments are in order. The brightness perturbations 
can be classified in terms of their transformation properties under rotations in the three- 
dimensional Euclidean sub-manifold. This means that all the brightness perturbations will 
have a scalar, a vector and a tensor contribution. In the framework of the ACDM scenario 
the tensor and the vector fluctuations of the geometry are totally absent. Therefore, since we 
want to compute the circular polarization in the minimal situation, we will stick to the scalar 
modes of the geometry which are connected with the scalar modes of the brightness perturba- 
tions. As implied by the concordance model the background metric is taken to be conformally 
flat with signature mostly minus, i.e. g 



{r) r n flu where r]^ = diag(l, -1, -1, -1) is the 
Minkowski metric. In Eq. f )2.3ip . and ijj represent the scalar fluctuations of the metric in 
the conformally Newtonian gauge, i.e. <5 s goo — 2a 2 </> and S s gij = 2a 2 5 i j'ip. Always Eq. f)2.3ip . 
v b is related to the baryon velocity (see appendix [SJ i.e. the center of mass velocity of the 
electron-ion system. In Eqs. fl2.3ip - fl2.34p e' denotes, as usual, the differential optical depth, 

Le - - 8 

<7 7 e = -Trr-g. (2.35) 



7c 



a Q 



11 Note the presence of the Doppler term arising in the source term of the intensity. The collisionless part 
of the evolution of the intensity perturbations is well known. For the interested reader it is derived, within 
the present conventions, in Ref. |51] and within different gauge choices. 
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After integration over cp' the source terms appearing in Eqs. fl2.3ip - fl2.34p can be written as 

3 n 



dfi' 



+ 
+ 2 



167T J-1 
P\\ ~ P\2 + Pi! — P22 



P n + P12 + P21 + P22 Ai(/i, n') + 
A Q (//,//) + 2 \p 13 + P 



23 



Pu + -P24 
3 n 



dfi' 



+ 
+ 2 



167T J-1 

Pll — P\2 — P2I + -^22 



A v (fi,fi' 

Pn + -P12 — P21 — P22 
A Q (/i,/i') 



(2.36) 



AiOx,//) + 



Cv(ji,u) 



Pl3 — P23 

3 f 1 

d/j,' 



Au + 2 



Pl4 — ^24 



A v (/i,/i') 



P31 + P. 



32 



167T J-1 

+ 2P 33 Au(/i, /i') + 2P 34 A v (/i, // 

3 rl 



Ai(/i,/i') + 



^31 " i 3 : 



32 



P41 + ^ 



12 



167T J-1 

+ P 43 Au(/i, //) + P 44 A v (/i, /i' 



Ai (//,//') + 



P41 - P 



12 



Aq (//,//) 



Aq (//,//') 



(2.37) 



(2.38) 



(2.39) 



where the generic matrix element appearing in Eqs. f l2.36p -f l2.39p is the integral over y?' of 
the corresponding element of the scattering matrix of Eq. (12. 3p . i.e. 



r2n ^2 

PijQi, fl!) = -r Sij(fl,IJ,',(p') dif'. 

Jo ri 



(2.40) 



For immediate convenience it is appropriate to write down the explicit form of the various 
matrix elements appearing in Eqs. fl2.36p - fl2.39p : 



Pu 




// 




= 7r{2A3( W )(l-A* 2 )(l- 


- /i' 2 ) + C 


2 (cW 2 [a 2 (u,) 


P12 




// 


n') 




AlM + /e( 


w)A|(w) 




Pl3 




// 


fjf) 


= 0, 








Pu 




// 


//) 








P21 




// 


//) 


= 7rC 2 M/i' 2 [A 2 (u;) + / c 2 


; W )A 2 M], 


P22 




// 




= 7rC 2 [A 2 M + / 2 (u,)A 2 ( 


w)], 




P23 




// 


//) 


= 0, 








P~2A 




[L 


//) 


= 2vr/ e HC 2 (u;)A 1 (u;)A 2 






P31 




// 


//) 


= P 3 2(w,//,/x') = P 3 3(W 


//, //) = 


P 34 (w,/i,/i') = 0, 


P41 




/' 


//) 


= 2vr/ c (u;)C 2 (u;)A 1 (u;)A 2 


M/i/i' 2 , 
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P 43 (u;, //,//) 



27r/ c (a;)C 2 (a;)A 2 (a;)A 1 (a;)/i, 
0, 



2vrC 2 (o;) 



A 2 ( W ) + /»A 2 (u,) 



/i/i 



(2.41) 



Inserting the results of Eq. (12.41 p inside Eqs. ( 12 . 3 6 1) — ( |2 . 3 9 [) the explicit expressions of the 
various source terms can be obtained. The details of this standard manipulation are reported 
in appendix IB1 (see, in particular, Eqs. (IB . 1 1) f |BT2]) ) . The final result is: 

1 



4 



A 



10 



2A 3 (o;)(l - + 2C 2 (u) Ai(u) + /»Aftu;) (1 + 



+ 



2As(u)(l - /i 2 ) - C» [Ai(u) + /»A 2 » ) (1 + ^) 

6^/ e (a;)C 2 (a;)A 1 (a;)A 2 (a;)(l + yU 2 )Avi 
1 



(2.42) 



A 



10 



2(1 - ^)A 3 M - 2C» A((u) + f \u)Ai(u) (1 - 



+ 



2Aa(w)(l - /i 2 ) + CV)(a?(w) + /»A 2 (u,))(l - /i 2 ) 
6z/ c HC 2 (a;)A 2 (a;)A 1 (a;)( yU 2 - 1)A V1 
C 2 MPi(/i) 



5,: 



/ e (o;)A2(a;)A 1 (a;) 



2Aiq — S$ 



3 

-?' 

2 



A 2 (u;) + / e 2 HA 2 V) 



A 



VI 



(2.43) 



(2.44) 



while Cxj(cu, fi) vanishes identically. In Eqs. (I2.42p . ( I2.43[) and (I2.44p the following important 
combination has been introduced, namely: 



Sy 



A12 + A Q2 + A 



(2.45) 



It must be remarked that Sp is the standard source term arising in the treatment of CMB 
polarization when the magnetic field contribution is ignored in the scattering process (see, 
e.g. (331 E21 [53] and references therein). The result expressed by Eqs. (I2.42p . (I2.43P and 
(I2.44p does hold to lowest order in the guiding centre approximation (see Eq. (lA.lip and 
discussion therein). As it will be discussed in HI the results derived so far improve the 
accuracy of the radiative transfer equations in the case when large-scale magnetic fields are 
consistently included in the discussion. 

The Faraday effect of the CMB is just a rotation of the linear polarization of the CMB 
and it does not involve the generation of any circular polarization. Faraday effect has been 
recently treated in greater detail by including various effects which have been neglected in 
the past (12] (see also [54J). In the case of the Faraday effect, first the linear polarization 
is generated because of the quadrupole in the intensity of the radiation field and then the 
polarization is rotated. Also Faraday rotation is treated often in the uniform field approxi- 
mation but without the explicit contribution of the magnetic field intensity to the scattering. 
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The present formulation improves also on the treatment of Faraday effect of the CMB (see 
[TO] for an introduction) even if, to keep the discussion self-contained, the focus will be on 
the generation of the V-mode polarization. 

In connection with the Faraday rotation it is appropriate to mention the different physi- 
cal nature of the approximations often employed in the discussion of large-scale magnetism. 
In the present paper, as already mentioned, the guiding centre approximation has been em- 
ployed. This approximation (see appendix |A|) is particularly sound in the case of scattering 
problems when the wavelength of the scattered photons is much shorter than the inhomo- 
geneity scale of the magnetic field [48]. The guiding centre approximation does not break 
explicitly the isotropy of the background. As already mentioned in the previous paragraph, 
Faraday rotation can be discussed in the uniform field approximation [53 ESI EH EH] ■ The 
uniform field approximation is independent upon the guiding centre approximation: indeed, 
for instance, in the studies of [55j ES E3 EH] the magnetic field does not contribute to 
the scattering matrix while it does rotate the polarization plane of the CMB. The uniform 
field approximation holds provided the magnetic field is not too strong. In the latter case 
a (new) preferred direction in the sky pops up; a potential correlation of the aj_ x m and 
aj +l m multipole coefficients is induced and, from this observation, uniform magnetic fields 
can be constrained [591 EQj EH E2]. This last case is not directly related to the present 
considerations. 



3 Circular polarization from magnetic knots 

It is appropriate to highlight a possible connection between the occurrence of circular polar- 
ization and the topological properties of the magnetic flux lines. The occurrence of circular 
polarization is directly related to the Lorentz force acting either on the individual charge 
carriers (i.e. v e> i x B) or on the Ohmic current (i.e. J x B). In a plasma characterized by 
a finite conductivity, the presence (or absence) of Lorentz force term can be directly related 
to the topology of the magnetic flux lines. The topology of the magnetic flux lines can be 
classified in terms of the so-called magnetic helicity, i.e. 

N B = I d 3 xA-B, (3.1) 

where A is the vector potential; Eq. (13. ip is the magnetohydrodynamical analog of the 
kinetic helicity, i.e. 

J\f v = J d?x v ■ uj, u = V x v (3.2) 

where uj is the vorticity. In simply connected domains, the magnetic helicity is gauge- 
invariant provided the normal component of B vanishes at the boundary surface of the 
integration volume V. Furthermore, the helicity is also gauge-invariant if the integration 
volume is given by a single (or multiple) magnetic flux tube. An important property is that 
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the magnetic helicity is conserved in a highly conducting plasma. In particular, it can be 
shown, that (see, e.g. [9] or the seminal paper of |47J) 



dH R 1 



dr 47TO" Jv 



d 3 xB • V x B, (3.3) 



where a is the conductivity. In the limit a — > oo the magnetic helicity is exactly conserved. 
In minimizing the total magnetic energy with the constraint that the magnetic helicity be 
conserved we are naturally led to the variational problem 



5 



d 3 x{ \B\ 2 -aA-B 

v 



0, (3.4) 



where a is a Lagrange multiplier. Since B = V x A, the variational problem of Eq. (13. 4p is 
equivalent to 

6 f v <Px(\V x A\ 2 -al- V x A^j =0. (3.5) 
By making the variation explicit, we have that Eq. (13.51) implies that 

VxVxl=oVxl (3.6) 
Going back to the magnetic field we have that the configurations 

VxB = aB (3.7) 

correspond to the lowest state of magnetic field energy which a closed system may attain. 
The variational approach leading to the condition (13 .7p is due to Woltjer [17] (see also, in this 
connection [131111] and [IUBE]). The configurations obeying Eq. (13.71) are closely related 
to the concept of magnetic knot [66] and may arise a consequence of the dynamics of the 
electroweak phase transition. The presence of pseudo-scalar interactions at the electroweak 
time can twist the magnetic flux lines of the hypermagnetic field and produce a primordial 
background of hypermagnetic knots [671 EE] • It has been speculated that baryogenesis can 
be related to the presence of hypermagnetic knots and a similar way of thinking has been 
pursued in [511 ED] • The themes discussed in [661 EZ] stimulated various investigations both 
on the dynamics of particles in hypermagnetic knot configurations [711 E21 E3] as well as on 
related ideas [7U [75] . 

In solar spectropolarimetry magnetic knots appear since the late sixties [76] (see also, for 
instance, [77J). An example of magnetic knot configuration is given by: 

B(z) = -B (sin kzx + cos kzy), (3.8) 

satisfying the force free condition discussed before and sometimes analyzed in connection 
with the polarization properties of the synchrotron emission [TBI EH]- The idea is to distin- 
guish the topology of the magnetic flux lines from the analysis of the circular polarization. 
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In other words: measurements of circular polarization can be used to infer not only the 
potential existence of large-scale magnetic fields but also their topological structure. As 
before we are considering here the situation where the ingoing Stokes parameters have no 
azimuthal dependence. The problem will now be to compute the various entries of the tensor 
Pij(fj,, fj/,z) which also depend upon the inhomogeneity scale of the knot. 

Consider, for simplicity, the case when protons are neglected and only the leading terms 
are kept in f e (uj). In this case we have, quite simply, that the phase matrix, after integration 
over ip', will be 

Pn(/x,//,z) = 7r[2(l-/i 2 )(l-/i ,2 )+/iV 2 ] 

+ TT/eH [3(1 - fl 2 )(l - /l' 2 ) + 1] COS [2(kz + </?)], 





z ) 


= 7T/i 2 + 7r{l+ fl 2 COS [2(kz + (f)]}f* (CO), 


P 13 (fi,fj,' 


z ) 


= 0, 


Pia{h, H 


z ) 


= 2irfifi'\Jl — fi 2 f c (u) sin [kz + <p\, 


Pi\{^1^ 


z ) 


= 7r/U ' 2 + 7r{l-cos[2(A;z + ^)]}/ 2 H, 




z ) 


= Tr + vrll-cos^-y?)]}/ 2 ^), 




z ) 


= 0, 


P 2 4(^,/i' 


z ) 


= 


Psiil^,^' 


z ) 


= -2n fif 2 {u) sin [2{kz + p)\, 




z ) 


= 0, 


Pw(n,v>' 


z ) 


= 0, 


P 34 (/i,/i' 


z ) 


= 2lX\j\ — j2 2 jJ,' cos [k z + ip]f e (u), 


P 4 l(/i,/i' 


z ) 


: -27^1 - fi 2 (fi' 2 - 2)/ e (w) sin [kz + <p], 




z ) 


= 2ix\j\- jj> 2 f c (uj) sin [kz + ip] , 


^43(^5 H 


z ) 


= 


-P44 (//,//' 


z ) 


= 2n^'[l + f 2 (u)]. 



By integrating over z the relevant matrix elements, i.e. 

r2Tz/k r2ir/k r2ir/k 

/ P 41 (n,n',z)dz = 0, / P 42 (n,fj,',z)dz = 0, / P 43 (fi,fi',z)dz = 0. 
Jo Jo Jo 

(3.9) 

The same integration, but applied to the linear polarizations, leads to a non- vanishing result. 
Indeed, the coupling of Aj to Aq is controlled by the following matrix element: 

Pll + Pl2 ~ P21 ~ P22 = 7T(1 - fi 2 ) ~ 37T/i' 2 (l - fi 2 ) 

+tt{3(1 - ^ 2 )(1 - //) + 2(1 + /1 2 ) cos [2{kz + p)]}f 2 (3.10) 
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whose integral over z does not vanish and is given by 



^(1 - /i 2 )(2 - 6// 2 ) + ^(1 - + 3// 2 + 2(1 - 3// 2 )]/ e 2 (3.11) 



The rationale for this occurrence stems from the fact that magnetic knots minimize the 
magnetic energy subject to the constraint the the helicity is constant. Indeed, over large 
scales, the minimization of the magnetic energy subjected to the constraint that the magnetic 
helicity is conserved is equivalent to the condition 

- - B 

Vx5=- (3.12) 

where L is has dimensions of a length and denotes the typical scale of the knot. Over very 
large-scales the displacement current can be neglected and, therefore, 

J — -?-V xB, J = en {vi - v c ) (3.13) 

47T 

But because of Eq. (13 . 1 2[) J = B/(47iL) and, therefore, J x B = 0. This shows that, rather 
generically, the vanishing of the Lorentz force implies, for these configurations, the vanishing 
of the circular polarization. 

There is another (indirect) way of appreciating this point. It is well known that, at finite 
conductivity and finite electron density it is possible to construct solutions of the Maxwell 
equations whose Poynting vector exactly vanishes both in the low frequency and in the high- 
frequency limit, i.e. E x B = 0. These solutions are often dubbed helicity waves since they 
do not carry momentum but rather helicity [80, ED E2] (see also [H31 El])- Consider first 
helicity waves in vacuo. A consistent solution of Maxwell's equations can be written, in this 
case, as: 

k 

E(z,t) = Bq — [sin (kz) sin (uor)x + cos (kz) cos (ur)y\, 

CO 

B(z,t) = £>o[sin (kz) cos (ur)i + cos (kz) cos (ur)y]. 

The above example can also be written as the superposition of circularly polarized waves 
propagating in opposite directions: 

-> k B 

E(z, t) = — — {[cos (kz — lot) — cos (kz + lvt)]x + [sin (kz + lot) — sin (kz — Lor))y}. (3.14) 
2to 

This class of solution can also be obtained at finite electron density and the pertinent dis- 
persion relations are 

4 = 1" 7 ^ ^ (3-15) 

where T is the collision frequency (i.e. the rate of interactions). The Ohm law can be written 
as: 

J = 7p ~^^vE, ^ = 7^- 3.16 

m e (1 — too) AttL 
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In the limit u — > the electric fields are suppressed by the conductivity while the magnetic 
fields will tend towards the force-free configuration recalled before. 

The configuration discussed here has some realistic features and the most relevant draw- 
back is that it is not localized in space. Localized knot configurations can however be con- 
structed (see [BB] and references therein and also [85]). It will be interesting to understand 
the scattering of photons also in these more realistic cases. In spite of that the physical mes- 
sage of the present exercise seems to be that circular polarization of the outgoing radiation 
is generic provided the underlying magnetic field does affect charged particles. As we saw 
such an inference is not automatic as long as knotted configurations maximize helicity but 
minimize the Lorentz force. In the latter case the scattering matrix might not be affected 
by the magnetic field if the correlation scale of the magnetic knot is much shorter than the 
Hubble radius at recombination. 



4 Estimates of the circular polarization 

Building up on the results of section [3] and taking into account the consideration of section 
H] it is now appropriate to solve the evolution equations of the brightness perturbations and 
to obtain explicit estimates of the V-mode power spectra. Since m e /m p <C 1 and f e (u>) <C 1, 
Eqs. (12.421) . (12.431) and (12.441) can be safely expanded in powers of (m e /m p ) as well as 
in powers of / e (w). The expansion of the scattering matrix in powers of (m e /m p ) is rather 
common (already in the absence of any magnetic fields) since, to leading order, the mean free 
path of the photons is chiefly determined by the scattering on the electrons. The expansion 
in (m e /m p )) is common practice in Boltzmann solvers (see, e.g. [85]). In the present case 
the same strategy will be employed by adding, however, a further expansion parameter, i.e. 

/.(")■ 

While the evolution equations of the brightness perturbations for the intensity and for the 
linear polarization have the first relevant correction going as fe{oo), the evolution equation 
for Ay has a source term proportional to / c (w). The three functionals appearing in Eqs. 
(I2.42p . (I2.43P and (I2.44p can then be expanded in powers of (m e /m p ) and f e (oj) with the 
result that 



Ci(u,/i) 



Cq(w,/x) 



A I0 1 + 



P 2 (» + 2 



2 + P 2 (/x) 



Avi + O 



^(/i) 



Sv 



^ 2 (/i) + 



S P - 2A M + 2iA V i 



O 



2 + P 2 (m) 



^{2/ e (a;)[2A I o - S P ) - + /»]A V1 



(4.1) 



(4.2) 
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(4.3) 



As anticipated, the source terms for the intensity and for the linear polarization have the 
first correction going as f 2 {uj) while the source term for the circular polarization starts with 
f e {oj). Higher order corrections to Eqs. (|4.ip . (14. 2[) and (14. 3[) can be computed, if needed 
recalling the results of Eqs. (I2.23[) — ( |2.25[) and of Eq. ( 12. 301) . Bearing in mind the results 
of Eqs. (14. ID . (14. 2 p and (14. 3p . to lowest order both in (m e /m p ) and in f e {u) the following 
system of brightness perturbations can be obtained: 



A[ + n^(Ai + <j>) + e'Aj = + e' 



/if b + Aio — S P 



A' P + n%A P + e'Ap = 
A^ + n%A v + e 'A v 



3(1 -fi 2 



Sv 



2A I0 — Sp 



3 A 

-zA Y1 



(4.4) 

(4.5) 
(4.6) 



where all the corrections 0(f 2 ) have been neglected. Note that in Eq. (14.51) Aq has been 
replaced by Ap, i.e. the brightness perturbation for the polarization degree P = V Q 2 + U 2 . 
In equivalent terms, as customarily done, we could have chosen the frame where Au = 0. 
Following the same notation, the source term Sp of Eq. (12.451) will become Sp = (Ai 2 + 
A P0 + A P2 ). If / e (w) = in Eq. (14. 6 p the standard set of brightness perturbations is quickly 
recovered. In this case the procedure will be to integrate the equations by assuming, for 
sufficiently early times, that the baryons are tightly coupled with the electrons implying 
that the baryon velocity is effectively equal to the dipole of the intensity, i.e. t>b — — 3iAn. 
This is, in a nutshell, the lowest order in the tight coupling expansion. To lowest order in 
the tight-coupling expansion the CMB is not polarized in the baryon rest frame, i.e. Ai 7^ 
An 7^ but Ai 2 = Ap 2 = Apo = 0. To first order in the tight coupling expansion (linear) 
polarization is generated and it is proportional, as expected, to the photon quadrupole which 
can be computed from the lowest order dipole. To summarize the approximations exploited 
so far we have that: 



• the scattering matrix has been derived in the guiding centre approximation; 

• the brightness perturbations have been then expanded for f e {oo) < 1; 

• the tight-coupling approximation is not invalidated by the new form of the evolution 
of the brightness perturbations. 

Before giving the details on the line of sight solution of Eqs. (I4.4p . (" 14 . 5 [) and (14. 6[) it is 
appropriate to pause a moment on the explicit numerical value of / e (^) 

^H^SX^X^!). (4.T) 
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For z ~ z ICC ~ 1091 (see e.g. / e (w) is of the order of 1(T 9 for nG field strength^. 

In Eq. (14. 7p 5 U denotes the uniform component of the magnetic field, i.e. we are assuming 
that the magnetic field is uniform since this is the simplest approximation in which that 
heat transfer equations can be analyzed. The considerations reported here in the uniform 
field approximation can be generalized to the case when the magnetic field is characterized 
by a given power spectrum. As already mentioned at the end of section [2] the uniform field 
approximation is more accurate in the present case than in the case of Faraday rotation 
which will be left for future discussions. 

The estimate of Eq. ( 14. 7ft can be further reduced by going to higher angular frequencies 
where, typically, nearly all CMB experiments are operating 13 !. Even if f e (co) can be rather 
small the question remains on the relative magnitude of the VV, VT and BB correlations and 
this will be one of the points discussed hereunder first for large angular scales (i.e. £ < 100) 
and then for smaller angular scales when dissipative effects are important. 



4.1 Line of sight solutions 



Equations (I4.4p . (14. 5p and ( 14. 6 j) can be solved, formally, by integration along the line of sight 
and the result of this step can be written, in Fourier space, as 



Ai(k,n,T ) 



+ 



K{r) 



Aio + fivb 



Sv 



3(1 



where, as usual, x 



Ap(fc,/!,To) 
Ay(k,H,U),T Q ) 

k(r — r) and 



4 



/C(r)5 P e-^ dr, 



2A I0 — S-p 



^Avi \ e-^ x dr, 



e(r, tq) = / x e n e adr , ^C(t) 



e e 



(4.8) 
(4.9) 
(4.10) 

(4.11) 



For large scales the visibility function /C(r) can be taken as sharply peaked at the recom- 
bination time. For smaller angular scales the (approximately Gaussian) width is essential 
to obtain sound semi-analytical estimates. Recalling the specific form of the lowest order 
Legendre polynomials [88l [89] 



PM 



Pi(/i) = n, 



(4.12) 



12 The dependence upon the redshift com.es about since the electron and proton masses break the Weyl 
invariance of the system (see appendix [K\ and, in particular, Eqs. (IA.4I) (| A.5|) . 

13 Just to have an idea the Planck explorer satellite is observing the microwave sky in nine frequency 
channels: three frequency channels (i.e. v = 30, 44, 70 GHz) belong to the low frequency instrument (LFI); 
six channels (i.e. v = 100, 143, 217, 353, 545, 857 GHz) belong to the high frequency instrument (HFI). The 
five frequency channels of the WMAP experiment are centered at 23, 33, 41, 61 and 94 in units of GHz. 
Neither WMAP nor Planck are sensitive to the circular polarizations. 
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Eqs. (|4.4p -( p~5|) and f)4.6p can be reduced to a hierarchy of coupled evolution equations for 
the various multipoles. Multiplying Eqs. f)4.4p . (14. 5[) and (14. 6 p by Po(^) = 1 and integrating 
over fj, between —1 and 1, the following relations can be obtained 



A; o + fcA n = ip', 
e' 

A' P0 + A;A P1 = -[A P2 + A I2 - A P0 ] 



A^o + kA 



VI 



i'A 



vo- 



(4.13) 
(4.14) 
(4.15) 



If Eqs. (I4.4l) - (l4.6p are multiplied by -Pi(/i), both at right and left-hand sides, the integration 
over fi of the various terms implies: 



-a; 



ii 



2 



5 A '» 



An + 7r:V h 



2 k 
-A' P1 - -kA P2 + -Ap 

2 k 
A' yi - -kA Y2 + -A vo = 



e'A P1 
3 



-A Y1 + -f c [2A l0 



A 



12 



A 



P2 



A 



PO 



The same procedure, using P2 (/■«), leads to: 
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-A; 2 - -fcAi 3 + -kA 

5 5 

3 2 
-A P2 - -A;A P3 + -fcApi 



9 1 

— A K (A P0 + A P2 ) 

10 10 V ' 



9 



A 



A V2 + 



A; 



3A V2 - 2A 



VI 



.10 
-e'A V 2. 



P2 



1 

To 



(A 



PO 



A 
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(4.16) 
(4.17) 
(4.18) 

(4.19) 
(4.20) 
(4.21) 



For i > 3 the hierarchy of the brightness can be determined in general terms by using the 
recurrence relation for the Legendre polynomials (see, e.g. [SSI EH]); the result for £ > 3 is: 

k 



A' u + t'A u = 
A' pe + e'A Pl 
A' yi + e'A ve 



2£ + l 
k 



- (^ + l)A I(m) ], 



o , fI [£A PM) -(l + l)Ap (w) ], 



(4.22) 



We are now ready to compute the evolution of the various terms to a given order in the 
tight-coupling expansion parameter r c = |l/e'|. After expanding the various moments of the 



brightness and of the baryon velocity in powers of r c , 

Ay£ = Ay^ + T c 5ye, 
Vb = V h + r c 5 Vh , 



An = A u + t c 8u, 
A P £ = A Pe + r c 5pe, 



(4.23) 
(4.24) 



the obtained expressions can be inserted into Eqs. fl4.13p - fl4.18p and the evolution of the 
various moments of the brightness function can be found order by order in r c . To zeroth 
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order in the tight-coupling approximation, t> 6 = — 3iA I1; while Eqs. (I4.14p and (14.1 7p lead, 
respectively, to 

Ap = A I2 + A P2 , A P1 = 0, A vo = 0. (4.25) 
Finally Eqs. (Q5]l and pD| imply 

g 

9Ai2 = Aq + Aq 2 , 9A Q2 = A qo + A I2 , A V i = -i/e(w)Ai . (4.26) 

Taking together the four conditions expressed by Eqs. (I4.25P and (I4.26P we have, to zeroth 
order in the tight-coupling approximation: 

A Qe = 0, £>0, A u = 0, £>2, (4.27) 
A w = 0, e^l (4.28) 

Hence, to zeroth order in the tight-coupling, the relevant equations are 

v b = -3iA n , (4.29) 
A~; o + kA u = if/, (4.30) 

A V i = ^/eMA I0 . (4.31) 

Even if to zeroth order in the tight coupling expansion we have that the linear polarization 
is absent. To get a non- vanishing linear polarization we have to go to first-order where the 
monopole and the dipole of the linear polarization are proportional to the quadrupole of the 
intensity; at the same order in the perturbative expansion the a non-vanishing quadrupole 
of the circular polarization is also generated. Indeed, recalling Eqs. (14 . 2 3 [) and (I4.24p . the 
first-order results can be written as 

5 1 

^QO = |^I2, $Q2 = -^>l2i (4.32) 
fa = jjrkKn, 5 V2 = j^Avi. (4.33) 
From Eqs. (I4.32p and (I4.33P the line of sight solution of Eq. (I4.10p can be written as 

A(k,fi,u,T ) = £° }C(r)^f e (oj)A lo (k,r) - ^(fcr c )A n }e-^dr, (4.34) 

where Eqs. (I4.29P and (I4.30p have been also used. The result of Eq. (I4.34p has been used in 
[87] and the present discussion corroborate and extends those results. 

In the sudden decoupling approximation the visibility function becomes effectively a 
Dirac delta function and the second term (proportional to the dipole of the intensity) can 
be neglected in comparison with the monopole whose evolution equation can be written as 

a "° + 74TT a :° + 3(74TI) A ^ e * + e - (435) 
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where = (3/4)p b /p 7 is the baryon to photon ratio and where 




(4.36) 



(4.37) 



Equation (I4.37P accounts for the presence of an inhomogenous magnetic field (see [Ml ESI [96] 
for a definition) but in the estimates we are going to present here, the inhomogeneities 
stemming from the magnetic field itself will be neglected for consistency with the uniform 
field approximation. 

It is finally appropriate to remark that it is possible to go improve on the accuracy of the 
line of sight solutions for the V-mode polarization. For this purpose the technique developed 
in Ref. [53] will be extended to the line of sight solution of the V-mode polarization. To be 
specific consider, again, the system of Eqs. ( 14. 8p . (14. 9 p and ( I4.1(jp . In particular we shall be 
interested in improving on Eq. (I4.34p which has been derived from Eq. ( I4.10p . Instead of 
using the lowest order tight-coupling solution for the linear polarization source it is possible 
to resum the perturbative expansion by solving an effective evolution equation for Sp. Recall, 
for this purpose, that the conformal time derivative of Sp, i.e. S P = (Aj 2 + A' P0 + A P2 ) can 
be expressed as the sum of the evolution equations of the separate multipoles. In particuylar 
A' P0 can be expressed from Eq. ( I4.14p . A' l2 can be expressed from Eq. (I4.19P and A' P2 from 
Eq. (I4.20p . Summing up the various contributions and rearranging the relevant terms an 
effective evolution equation for Sp can be obtained and it is 



which can be solved in different ways. For instance, if only the intensity dipole is kept, at 
the right hand side of Eq. (14.381) the result for Sp is 



If we now insert Eq. (14.391) into Eq. (14. 10p the integrals can be performed, up to some point, 
with analytic techniques. This result improves then Eq. (14 .34 p . 

The considerations developed so far suggest the following physical picture. To lowest 
order in the tight-coupling expansion the presence of a magnetic field produces a dipole of 
the circular polarization. The dipole of the V-mode can be fed back into the line of sight 
solution to obtain the higher multipoles in full analogy with what is customarily done in 
Boltzmann solvers. While the circular polarization is building up from the monopole of the 
intensity, the dipole of the intensity sources, to first-order in the tight coupling expansion, 
the linear polarization and, in particular, Apo and Api. The interesting aspect of this 
analysis is that, indeed, to lowest order in the tight coupling approximation the CMB is 
circularly polarized if a pre-decoupling magnetic field is present. The linear polarization 




(4.38) 




(4.39) 
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is generated to first-order in the tight-coupling expansion. At the level of the amplitude, 
as it will be shown, the angular power spectrum of the V-mode is always smaller than the 
E-mode spectrum which arises directly from Ap. The reason for this occurrence stems from 
the specific value of f e (u>). 



4.2 Large-scale limit 

For large angular scales the circular polarization will then be given by 

Av(M,To) = ~ r JC(r)f e (co)e-^- T ^A m (k,r)dT, (4.40) 
3 Jo 

where Aio can be determined from Eqs. (I4.29l) - fl4.3ip and (14.351) . The terms arising in Sp 
have been neglected since they vanish to lowest order in the tight coupling expansion. By 
assuming that /C(r) is a Dirac delta function centered at recombination (sudden decoupling 
approximation) we shall have that 

g 

A y {k,(i,To,u) = -f e (u)e~ ifiX ^A l0 (k,r Iec ). (4.41) 
Let us now compute the angular power spectrum of the circular polarization. 

A v (n,r ) = 5>2r, m (™), (4-42) 

thus we will also have that 

affi = J dhA y (h, r )Y; m (n) = J dh J d 3 kA Y (k, /a, r ) Y? m (h) (4.43) 

Recalling the explicit expression of Ay(k,fj,, r ) 

(V) _ 8jeH 



dh J rf 3 fcF; m (n)/ie-^A I0 (A;,r r( 



3(2tt) 3 / 2 

J d nJ d 3 kY; m (n)^(e-^A m (k } r tec ). (4.44) 

The integration over dh = smddipdd = —d^dtp can be performed in explicit terms since 
e -fyiB can expanded in Rayleigh series and the final result is 

«S = -^-o^p/ e V4Tv / 27TT / d*k(^jA 10 (k,T Tec ). (4.45) 
The angular power spectrum of the circular polarization can then be written as 

rn=—£ 

2567T „ 9/ , r°° dk k 3 fdji\ 2 . A ., ..„ 
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where the dependence on the angular frequency u has been explicitly included in the ex- 
pression of the angular power spectrum. The large-scale estimate of Ai (k, r rec ) follows from 
Eq. (I4.13P by neglecting the dipole which is negligible for those wavelengths which are still 
larger than the Hubble radius around the redshift of recombination: 

AioO,T rcc ) = A I0 (A;,t*) +ip(k,T Tec ) -ip*(k), (4.47) 

where, by definition, ip(k,T*) = ip*(k) and 0(fc,r*) = 4>*{k) are the values of the metric 
fluctuations at r* <C r eq . For r ~ kr* <C 1 and, in the minimal ACDM scenario, the initial 
conditions are (predominantly) adiabatic, i.e. 

,{k) 
2 

implying 

Ai (fc,r rcc ) = ^±^ 7^(fc), Mk) = (l + \Ru)Uk), (4.49) 

where FL is the fractional contribution of the massless neutrinos to the radiation back- 
ground_|. The angular power spectrum is 

r2 / r? j_ i k \ 2 / u\ n.,-1 



6 7 (k, n) = -2^(fc), i/}(k, n) = ip*(k) K*{k) = -if>.(k) - ^tt 1 , (4.48) 



(W), v 512tH / ^+15 \ 2 2 /fcoX"- 1 {vv) 

da? 



it Jo \( 



The derivative of the spherical Bessel function can be expressed in terms of the appropriate 
recurrence relations, namely [88| 189]: 



-f = -3t(x) - j t+1 (x). (4.50) 
ax x 



The integral becomes then: 
X^\£ } n s ) = £(£ + I) / da; a; 



„«s— 5 



^ Jf+1/2{ X ) + x2 ^£+3/2( X ) — 2£xj£ +1 / 2 (x) Jl+3/2(%) , 

(4.51) 

which can be explicitly computed [88, 89J. The final result for the angular power spectrum 
can therefore be written as 

£(£+l)„ (yY)/ v 256vr/ R u + 15 \ 2 „ _ /^vi 



2tt 

where 



- iH^J '»MS) " l(vv, <«- <« 2 » 



£(£ + 1)[4£(£ + 1) - (n s - 1)K - 2)(n s - 4)]T(3 - n s )rf£ - f + f 
jf V) K) = 7 ^- ^- ^-A (4.53) 

2 6-n sr / 2 - f JW 3 - f Jrf | - f J 



14 Neutrinos are taken to be massless for consistency with ACDM paradigm even if the effect of the masses 
could be included without appreciable changes in the forthcoming numerical estimates. 
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Following the same technique we can estimate the cross-correlation between temperature 
and polarization as: 



^±^ v » = — /cM ^ +15 )^- 15 v^y s ' 1 x^)(£, 

2tt 1 K ' 75 J K ' (4^ + 15)2 \k p J v ' 



where 



Jo 



£ J e+i/2( x ) - xJe+i/2(x)J e+3/2 (x) 



(4.54) 



(4.55) 



By integrating the above expression we have that: 



+ l)(2-n s )r(^2-f)r^+f - l) 



4^r |-f r 



f + 3 



(4.56) 



In Fig. [TJ the VV and the VT angular power spectra are illustrated for large angular scales 
(i.e. I < 100). Both in Fig. [Hand in Fig. [2] a double logarithmic scale has been used. In 
Fig. [1] the uniform magnetic field intensity is fixed while the frequency ranges between 10 
GHz and 30 GHz. Note that 30 GHz corresponds to the lower frequency band of the Planck 
explorer satellite which is unfortunately not sensitive to the circular polarization. In Fig. [1] 
the cosmological parameters have been fixed as 

(fi b , tt c ,tt de , h , n s , e rc ) = (0.0441, 0.214, 0.742, 0.719, 0.963, 0.087), (4.57) 

corresponding to the best fit of the WMAP data alone in the light of the vanilla ACDM. In 



WMAP 5-yr data alone, B = 1 nG 



WMAP 5-yr data alone, B = 1 nG 
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Figure 1: The VV and VT angular power spectra are illustrated for large angular scales and 
for different frequencies at a fixed value of the magnetic field intensity. In both plots a double 
logarithmic scale has been employed, i.e. on both axes we plot the common logarithms of 
the corresponding quantity. 
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WMAP 5-yr data alone, v =10 GHz WMAP 5-yr data alone, v = 1 GHz 
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Figure 2: The VV and VT angular power spectra are illustrated for large angular scales and 
for different values of the magnetic field intensity at a fixed value of the frequency. As in 
Fig. [1] a double logarithmic scale has been employed in both plots of the present figure. 

Fig. [2]the frequency is fixed to 10 GHz while the magnetic field intensity changes from 1 to 10 
nG. The values of the magnetic field intensity are motivated by a recent analysis [34] where 
using the TE and TT correlations the parameters of a putative magnetized background are 
scrutinized. According to [33j El] nG magnetic fields and blue magnetic spectral indices are 
fully compatible with the measured values of the TT and TE angular power spectra. 

The results of Fig. [1] and of Fig. [2] show already that the values of the VT correlations 
are close to the magnitude of the BB angular power spectra from gravitational lensing as 
well as to the BB angular power spectra expected from the tensor modes of the geometry. 
It is important to remind here that the VV and VT angular power spectra are a direct 
consequence of, both, the weakly magnetized plasma and the adiabatic mode of curvature 
perturbations. Absent one of these two essential ingredients the net result would vanish. 
The considerations reported so far complement some of the results already reported in [87J. 
The large-scale (analytical) estimates will now be corroborated by the numerical results at 
smaller angular scales. 

4.3 Small-scale limit 

The visibility function vanishes for r 3> T rec and has a maximum around recombination, i.e. 
when 

e" + e' 2 = 0, ^ = -x c (r)n c (r)a Th a(T) = -e'. (4.58) 
ar 

The second expression in Eq. (14.581) clarifies that a minus sign appears in the time derivative 
of e(r, To) since r appears in the lower limit of integration. When the finite thickness effects of 
the last scattering surface are taken into account the visibility function can be approximated 
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by a Gaussian profile centered at r rec , i.e. 



/C(r) 
A/"(a rcc ) 
erf (z) 



Af(a rec )e 
r 2 1 



(T-Trec) 



/C(r)dr 



7T (7 

— / 



erf 



T~0 - T Ti 



y/2, 



a r 



+ erf 



(Tree \ 
y/2a ICC J . 



~ l dt. 



(4.59) 
(4.60) 
(4.61) 



The overall normalization M{a rcc ) has been chosen in such a way that the integral of /C(r) 
is normalized to 1: the visibility function is nothing but the probability that a photon last 
scatters between r and r + dr. Equation (I4.60p simplifies when t 3> T rec and r 3> cr rec , since, 
in this limit, the error functions go to a constant and Af(a rec ) — » cr~,i J2/tt. In the latter 
limit, the thickness of the last scattering surface, i.e. a rec , is of the order of r rec . The Gaussian 
approximation for the visibility function has a long history (see, e.g. j9(Jl ED E2J |93]). The 
WMAP data suggest a thickness (in redshift space) Az rec ~ 195 ± 2 which would imply 
that <T rec , in units of the (comoving) angular diameter distance to recombination, can be 
estimated as cr rcc /r ~ 1.43 x 10~ 3 . By using the finite width of the visibility function we 
have that 

at ,-\M „ , J- x i- 2 

a (V) - -5 



3(2tt) 3 / 2 ' 



,/ e \/W2£ + 1 / d 3 k 



dji 



dx 



e k t ji A M (/i;,r rec )e" 



(4.62) 



where k t = (v6/cr rec ), x = A;(r — r roc ). The visibility function, more realistically, has also 
a second peak at the reionization epoch (i.e. for z re = 11 ± 1.4). Also in this case the 
visibility function can be approximated with a Gaussian profile centered, this time, around 
r re and this consideration introduces a suppression going as e~ Erc . In the limit tq ^> r rec and 
to 3> °rec the integral to be computed is, therefore, 



+ 1) 



2tt 

Recalling that 



a 



CO 



8/eV) 



2tt s 



\A l0 {k,r Tec )\ 2 (i+l)£ 



dh\ 2 
dx) 



fdot\ 2 = 




\dxj 


X 2 



1 d 2 



(4.63) 



(4.64) 



the Bessel functions can be estimated in the limit of very large multipoles. A standard 
calculation then leads to 



l(l + l) c (w) 
2vr 1 
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fc 



n s — 1 



W 



n s — 5 



\/w 2 — ldw x 



x 



2/4( W ,£)e- 2 ^ 2 "' 2 + M^w,^^- 2 ^' 2 '" 2 ) 

+ £)c sb cos [2( 7A ^)] e - 2( ^ s)2w2 

+ AL n (w4)M n (wJ)^cos( lA wt)e-W e ^ + W^ w2 



(4.65) 
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where c s b = [3(i?b + 1)] 1//2 is the baryon sound speed. Following Ref. [51] we introduced in 
Eq. (14651) the quantified 



L n (w, £) = a n - P n In [wq e ], M n (w, t) = ot n + fi n In [wq e 



200 d A (z 




6 ' 

ISr = — ln(96), ^ w = — . (4.66) 

In Eq. f)4.66p the quantity c^z) is related to the angular diameter distance D A (z) as 
cIa(z) = y/VlMoHoDA(z)/2. Furthermore 7a and £r> can be estimated as follow£§ 



dl\z) 
7a = , = m 

'3R h (z + 1) 



(B = kD Q A{z)= (^) 5/4 W. (4.67) 



The Silk multipole is just expressed in terms of £ t and £^> as £s = ^T>/y^t +^r>- With 
the same approach it is possible to express, for small angular scales, the cross-correlation 

For more reliable estimates at small scales fully numerical methods should be employed 
and the results, consistent with the previous analytic estimates, are illustrated in Fig. [3jand 
IU In Fig. [3jthe VV and VT angular power spectra are illustrated on a double logarithmic 
scale. The spikes appearing in the VT correlation are the usual feature displayed when 
plotting the modulus of the cross-correlation. The same spikes occur when plotting the 
logarithm of the modulus of the TE power spectrum. In Fig. [3jthe magnetic field is fixed 
to 1 nG while the frequency of the channel is allowed to vary. In Fig. H] the VV and VT 
correlations are illustrated for a fixed value of the frequency but for various magnetic field 
intensities. 

The moment has now come to compare the signal from circular dichroism with the 
signals expected from the linear polarization. The spirit of the forthcoming considerations 
is just to compare the orders of magnitude of the different contributions. While this step 
is mandatory within the present analysis it is also not conclusive, from the experimental 
point of view. Indeed different signals and different correlation functions inherit from nature 
different systematic effects. The latter problem is of course extremely important and will 
not be treated here. 



15 It is understood that all the quantities are computed for at z — z rcc ; furthermore = pn/pM = 
4.15 x 1CT 2 wm(.z/1000). 

16 Following the usual convention we shall denote ujx = h^ilxo for a generic species. 
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WMAP 5-yr data alone, B = 1 nG WMAP 5-yr data alone, B = 1 nG 
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Figure 3: The VV and VT angular power spectra are illustrated for small angular scales and 
for different values of the frequency and for fixed magnetic field intensity. As in Fig. [2] a 
double logarithmic scale has been employed in both plots. 



WMAP 5-yr data alone, v = 10 GHz WMAP 5-yr data alone, v = 10 GHz 
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Figure 4: The VV and VT angular power spectra are illustrated for small angular scales and 
for different values of the magnetic field at a fixed value of the frequency. As in the previous 
figures, a double logarithmic scale has been employed. 

Let us first of all compare the V-mode signal with the measurements (and expectations) 
of the ACDM paradigm. A simple comparison is illustrated in Tab. [T] where the different 
angular power spectra are reported at the peak. The TT power spectrum at the first acoustic 
peak is of the order of 6000 (/iK) 2 . The TE and the EE angular power spectra peak for 
larger multipoles. The WMAP team measured with reasonable accuracy the region of the 
first anticorrelation peak of the TE power spectrum and the lowest multipoles of the EE 
power spectrum. The recent QUAD measurements gave a rather interesting evidence of 
the oscillations in the EE angular power spectra at larger £. The absolute value of the TE 
angular power spectrum peaks around I ~ 750 and it is of the order of 130 (fiK) 2 . The 
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EE correlation reaches a value of 40 (/iK) 2 for i ~ 1000. The quoted figures are consistent 
with the expectations of the ACDM scenario with no tensors (also sometimes called vanilla 
ACDM). In the vanilla ACDM the only potential source of B-mode polarization is represented 



Data 


^peak 


WW + 1)C£2/(2tt) 


TT 


220 


5756 (/xK) 2 


EE 


1000 


40 (/iK) 2 


TE 


750 


130 (/iK) 2 


VT 


i < 50 


10- 6 (/iK) 2 


VV 


t < 50 


10~ 14 (/iK) 2 



Table 1: The values of the different angular power spectra at the peak (llustrative figures 
for B u = InG and v = 10GHz. 

by gravitational lensing of the primary anisotropics. The typical values of the induced 
BB angular power spectrum range between 10 -8 (/iK) 2 and 10 -5 for £ < 50. In Tab. [1] 
the expectations for the VV and VT angular power spectra are reported in the case of a 
hypothetical low frequency instrument sensitive to V-mode polarization operating in a band 
with central frequency of the order of 10 GHz. The intensity of the (comoving) magnetic 
field has been taken to be 1 nG. As previously discussed both analytically and numerically 
the VV and VT power spectra are larger for low multipoles. For larger multipoles, however, 
the effects of the thickness of the visibility function and of the diffusive damping come into 
play only for I ~ 1000. In this sense the range of multipoles highlighted in Tab. [D should 
be complemented with the results illustrated in Figs. [H 121 El and SJ 

The comparison summarized in Tab. [1] does not contemplate the BB angular power 
spectrum stemming from the tensor modes of the geometry which is regarded as the main 
target of running experiments such as Planck. The B-mode power spectrum induced by relic 



Data 


r T 


n s 


n A 


^MO 


WMAP5 alone 
WMAP5 + Acbar 
WMAP5+ LSS + SN 
WMAP5+ CMB data 


< 0.43 

< 0.40 

< 0.20 

< 0.36 


0.986 ±0.22 
u.ao<j_ 020 
0.968 ±0.015 
0.979 ± 0.020 


77n+ - 033 

u - ' ' u -0.032 

0.767 ±0.032 
0.725 ±0.015 
0.775 ± 0.032 


n oonO.032 
U.ZOU_ 033 

0.233 ±0.032 
0.275 ±0.015 
0.225 ± 0.032 



Table 2: The change in determination of the parameters of the tensor background for three 
different choices of cosmological data sets. 

gravitons peaks for I ~ 90 corresponding to an angular scale of roughly 2 deg. The signal, 
however, depends upon rx (i.e. the ratio between the tensor and the scalar power spectrum) 
for which only upper limits exist. Depending upon the data sets chosen for the analysis the 
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putative limit on slightly changes. The situation is quickly summarized in Tab. [2] where 
the upper limits on r<i are reported at the pivot scale k p = 0.002 Mpc^ 1 and in the case 
where the scalar spectral index does not run. Slightly larger values of rx are allowed if n s 
is allowed to run but this aspect will not be essential for the present considerations The 



WMAP 5-yr data alone WMAP 5-yr data alone 
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Figure 5: The VV and VT angular power spectra are compared with other polarization 
signals arising in the ACDM paradigm and in its neighboring extensions. 

comparison between the VV and VT power spectra and the other polarization power spectra 
is also illustrated, more visually, in Fig. |5j In the plot at the left of Fig. [5] the upper curve 
corresponds, as indicated by the legend, to the EE angular power spectra obtained from the 
best fit to the WMAP 5-yr data alone (see Eq. ( 14.57P ). In the plot at the right the B-mode 
polarization is illustrated when a tensor mode contribution is allowed. 

It is tempting to speculate, at this point, that, indeed, low frequency instruments could 
make the difference for scrutinizing a potential V-mode polarization. In this respect the 
results and the techniques of [381 EH HO] (as well as the earlier results of [361 EZ]) could 
be probably revisited in the light of the considerations developed here. It has been shown 
that the VT correlation for a comoving magnetic field from 5 to 10 nG can be as large 
as 10~ 5 (/iK) 2 at 10 GHz for £ < 20 (i.e. large angular separations). This means that for 
frequencies (9(MHz), the resulting signal could be even 6 or 7 orders of magnitude larger 
than a putative B-mode signal from gravitational lensing which is between 10 -8 (/iK) 2 and 
10" 6 (/iK) 2 . 

5 Concluding considerations 

In the present paper it has been argued that the presence of a large-scale magnetic field prior 
to equality can affect the photon-electron and the photon-ion scattering. In this process the 

17 In the case of running the bounds on rr range from 0.58 (in the case of the WMAP 5-yr data alone) to 
0.54 if we combine the WMAP data with the large-scale structure data and with the supernova data. 
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radiation becomes circularly polarized and the induced VV and VT angular power spectra 
have been computed. The analysis reported in the present paper has to be regarded as very 
preliminary. Indeed the considerations reported here can be refined both at the theoretical 
as well as at the more observational level. At the same time this investigation certainly 
opens the way for a more direct use of the circular polarization as a specific diagnostic of 
pre-decoupling magnetism. The reported results might also be regarded as a modest spur 
for those observers and experimenters whose aim is a direct measurement (or a plausible 
upper limit) on the circular dichroism of the Cosmic Microwave Background. 
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A Derivation of the scattering matrix 



In a conformally flat background geometry characterized by a metric tensor g^v{r) = a 2 (r)?7 Atl/ 
the electron-ion plasma can be described by the well known set of two-fluid equations (see, 

e.g., mm)- 

V ■ E = 4ire(m - n e ), V ■ B = (A.l) 

OE OB 

V x B = Aireimvi - n e v e ) + — , Vx£ = -— , (A.2) 

or or 

where the electromagnetic fields as well as the concentrations of electrons and ions are 
comoving, i.e. 

B(x,t) = a 2 (T)B(x,r), E(x,t) = a 2 (T)£(x,r) n eti (x, r) = a 3 (r)n ej i(x, r). (A. 3) 
The evolution equations of the electron and ion velocities can be written as 

^ + nv e = -—\E + v e x B], (A A) 

ar m c a 

dv~ & —* — * 

-± + Hv i = [E + ViXB], (A.5) 

ar m p a 

where m e and m p are, respectively, the electron and the ion masses; the velocities are related, 



as usual, to the comoving three-momentum v^i = q e ,i/ y <?e,i + m e, P a2 - The explicit depen- 
dence upon the scale factor at the right-hand side of Eqs. (IA.4[) and (IA.5[) arises because the 
plasma is cold: both electrons and ions are non-relativistic and the mass dependence breaks 
explicitly the Weyl invariance of the whole system already at the level of the Vlasov-Landau 
equations for the distribution function (OH IS] . Finally, in Eqs. ( 1A.4I) and (1A.5I) H = a'/a 
enters directly the Friedmann-Lemaitre equations: 

H 2 = ^a 2 p u H 2 -U' = ArrGa 2 ( Pt + p t ), (A.6) 

where p t and p t denote the total energy density and the total pressure while the prime stands 
for a derivation with respect to the conformal time coordinate r. Since the electron and ion 
concentrations are comoving, they simply obey the following pair of equations 

< + V • {n e v e ) = 0, n[ + V • (71$) = 0. (A.7) 

Further details on the description of globally neutral plasma in Friedmann- Robertson- Walker 
backgrounds can be found, for instance, in [T0J [33l [3H [91] (see also [95], [96] for earlier results). 

The purpose is now to derive the components of the scattering matrix Sij(uj, /i, if, //, <p') 
introduced in Eq. (12.31) . In the dipole approximation [JJ, the outgoing electric field can be 
written as: 



E out {x,r) = — 



, . d.J 
r x r x — 
dr 



(A.8) 
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where 

J(x, t) = e(niVi - n e v e ) = en (vi - v e ); (A. 9) 

the quantity n = n e = n\ denotes the common value of the (comoving) electron and ion 
concentrations. In components the outgoing electric field can be written as 

E k = —[(J'-f)r k -J' k ], (A.10) 
nor 

where the prime denotes, as usual, a derivation with respect to the conformal time coordinate. 
The maximum of the microwave background arises today for typical photon energies of the 
order of 10 -3 eV corresponding to a typical wavelength of the mm. At the time of decoupling 
(i.e. z dcc — 1089 according to [221 E31 El]) the wavelength of the radiation was of the order 
of 10 _3 mm ~ /im. Since the magnetic field we are interested in is inhomogeneous on a much 
larger length scale we can use the guiding centre approximation |48j stipulating that 

Bi(x, r) ~ Bi(x , t) + (x j - xfodjBi + ... (A.ll) 

where the ellipses stand for the higher orders in the gradients leading, both, to curvature 
and drift corrections which will be neglected throughout. Fixing a local coordinate system 
with three orthogonal axes x, y and z the components of the accelerations for the electrons 
are 

^ = -^ °? = & E - (A - 12) 

where / e , i = c, i/w; wb e , i and w p e , i are, respectively, the Larmor and the plasma fre- 
quencies either of the electrons or of the ions (see Eqs. ( I2.27P and 12. 28]) ). The magnetic field 
is oriented along the z axis and only the lowest order in the gradient expansion is kept. Of 
course, as it is well known, higher order will induce both gradient drifts as well as curvature 
drifts (see, e.g. j3H]). We consider these terms to be negligible in the first approximation. 
Here we are interested in the scattering of electrons and photons not in the stellar atmosphere 
but rather at the decoupling time when the physical wavelength of the photons is minute in 
comparison with the inhomogeneity scale of the magnetic field which is of the order of the 
Hubble radius and even larger [331 El]- The guiding centre approximation is pretty robust 
as far as the magnetized scattering in concerned. For sufficiently small angular scales (i.e. 
£ ^> 100) the radial direction, in spherical coordinates, coincide (approximately) with the z 
direction. For this reason, in some related paper the modulus of the magnetic field has been 
taken &sB = z- B~h-B where the last approximate equality follows in the limit of small 
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angular scales. In the limit £ 3> 100 the two-sphere actually degenerate into a plane and this 
is the reason why one can trade the spherical decomposition for the plane wave expansion. 

It is also possible to write the difference of the accelerations of electrons and ions, namely, 
A = (a* 6 ) — a^) whose components in the local frame are 

,2 



A n 



A, 



A, 



P e 



47rn e[l-/e 2 M] 



P e 



47m e[l - P e (u)\ 



[A 1 {uj)E x -if e {uj)A 2 {uj)E y 
[Kx{u)E y + if e (u)A 2 (uj)E x 



U). 



P e 



A s (uj)E z 



(A.13) 
(A.14) 
(A.15) 



where Ai(w), A 2 (u) and A 3 (u) have been already defined in Eqs. (12T23]) . (12T241) and (12T251) . 
To compute the elements of the scattering matrix it is preferable to pass from the Cartesian 
components of the incident electric fields to the polar components; the relation between the 
fields in the two basis is given by: 



E x = E'fl cos i?' cos <p' — sin ip'E, 
E y = E'fl cos sin <p' + cos <p'E' t 
E 7 = -siw&'EL 



> 771/ 



(A.16) 
(A.17) 
(A.18) 



To avoid a proliferation of superscripts in the intermediate expressions the ingoing electric 
fields Efl and E' v are renamed as Ei and E 2 , i.e. 



Ei = EL E 2 = E' 



(A.19) 



Using Eq. (1A. 12j) the scattered electric field, in the dipole approximation, become 

E = Ej + E 9 0, (A.20) 
where recalling the notations of Eqs. (I2.23l) - fl2.26l) 



E#(u),ii,ij',<p,<p') 



V 



Ai(o;)C(w) /if [i'Ei cos ((f' — ip) — E 2 sin (cp' — cp) 



i-^Vi-AM^i 



iA 2 (u)f e (u)((u)ii 



jj'Ei sin (ip' — if) + E 2 cos (ip' — (p) j, (A. 21) 



C(^)AiH 



li'Ei sin {cp' — ip) + E 2 cos ((p' — ip) 



+ i/ c (w)A 2 (w)((w) 



ji'Ei cos {ip' — ip) — E 2 sin (ip' — ip) 



. (A.22) 



Thus, collecting the different factors and rearranging the final expressions we shall have that 
the scattered electric fields can be written as 



E 6 (u,fjL,f/,(p,(p') = -^[AEx + B E 2 ] 



E v (w,fi,fi',<p,<p') = -[CEt + DE^ (A.23) 
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where 

A(u,ii,ii',(p,(p') = CM/^'AiH cos (ip' - y?) - yl - - /i' 2 A 3 (w) 

- iA 2 {u)f e (u)C(u)fifi' sin (p' - (p) (A.24) 

B(u, [L, //, v?, V 2 ') = -((w)Mi(w) sin (<£>' - <p) - iA 2 (u)f e (u)C(u)fM cos <y - y>)(A.25) 

C(u,fi,fi',ip,(p') = /i'C(w)Ai(w) sin(y -</?) + i/ c (w)A 2 (w)C(w)/i / cos (</?' - y?),(A.26) 

D(u,/i,ii',(p,(p') = £(u)Ai(u) cos (y?' -93)- i/ c (w)A 2 (o;)C(w) sin (<p'-<p). (A. 27) 

The Stokes parameters of the scattered radiation can be related to the Stokes parameters of 
the incident radiation in terms of the appropriate scattering matrix S. As already mentioned 
in connection with Eq. (12.31) . arranging the outgoing and the ingoing Stokes parameters in 
a pair of column vectors 

X out = I v , U, V), X in = (Ji, J 2 , [/', ( A - 28 ) 

the outgoing Stokes parameters are given as X out = SZ in , where the various components of 
S are: 

5n = %^| 2 , <S 12 = %5| 2 , 

«S 13 = ^(A*B + BM), S 14 = % ^(A*B - B*A), (A.29) 
<S 21 = ^|Cf , <S 22 = -f |Z5| 2 , 

S23 = ^(C*D + D*C), S 24 = t^(C*D - D*C), (A.30) 

2 2 
<S 3 i = + AC*), S 32 = ^(B*D + D*B), 

S 33 = ^(A*D + AD* + BC* + B*C), 

■ 2 

S U = ^(A*D-AD* + BC*-B*C), (A.31) 
2 r 2 

• 2 • 2 

5 41 = l -^(AC* - A*C), 54 2 = ^(BD* - B*D), 

■ 2 

<S 43 = l -^-(AD* -A*D + BC* - B*C), 
2 r 2 

S 44 = ^-(A*D + AD* - B*C - BC*). (A.32) 
2 r 2 

Using Eqs. (IA.24p - flA.27p inside Eqs. flA.29p - flA.32p the various matrix elements can be 
readily obtained in explicit terms and have been reported from Eq. (12.70 to Eq. f)2.22p . 
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B Details on the derivation of the source terms 



It is appropriate to give few details on the derivation of the source terms reported in Eqs. 
(1^421) . ( EZ3D and (EE}. Using the results of Eq. ( 12H1 the source terms of Eqs. (12T3TD - 
f)2.34p can be rewritten in explicit terms as 



16 



2A 3 M(i-^ 2 )z 1 + c 2 M 



A?( W ) + / e V)A 2 (u,) 



(l + /^ 2 )2 2 



2A 3 (c)(l - /, 2 ) - C» (A 2 (c) + / c 2 MA 2 (u;)) (1 + ^ 2 ) 



+ 

3 



2, 



C Q ( W)A i) = ^{2A 3 (o;)(l- A i a )^-CV)(A?(a;) + ^(a;)A»( W ))(l-A* 2 )^ 



2A 3 ( W ) + ( 2 M U?M + / C 2 MA^) 



+ 



+ 4/ c HC 2 (a;)A 2 (a;)A 1 (u;)(/i 2 -l)Z 4 



(l-/i 2 )^3 



C v (a;,//) = ^|/// e C 2 (a;)A 2 (a;)A 1 (a;) 

+ /xC 2 M(A 2 (a;) + / c 2 (a;)A 2 (a;))z 4 }. 



2^5 — 2s 



(B.l) 



where the four quantities 2, with z = 1, 2, 3, 4 denote the integral of the various brightness 
perturbations, i.e. 



Z x = / 1 (l-/i ,2 )A I ^ / = ^[A I0 + A I2 ], 
J -i 3 

Z 2 = [ (1 + /i /2 ) Aj d//, = ^A I0 - ^A K , 
j-l 3 3 

Zs = Jjl - ^) A Q d/i' = ^[A Q0 + A Q2 ], 
2 4 = y // Ay g?// = — 2i Avi. 



(B.2) 



To derive Eq. HB .2j) the standard multipole expansion for the brightness perturbations has 
been assumed, i.e. 

A x (n, r) = £(-z)' (2£ + 1) A»P,(/x). (B.3) 



Using the explicit expressions of Eq. ( 1B.2I) inside Eq. ( IB. II) . the expressions reported in Eqs. 
(g35D and ( jXBj) are quickly recovered. 
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